Dynamics of karyotype evolution.

In the evolution of species, the karyotype changes with a timescale of tens to hundreds of thousand years. In the development of cancer, the karyotype often is modified in cancerous cells over the lifetime of an individual. Characterizing these changes and understanding the mechanisms leading to them has been of interest in a broad range of disciplines including evolution, cytogenetics, and cancer genetics. A central issue relates to the relative roles of random vs deterministic mechanisms in shaping the changes. Although it is possible that all changes result from random events followed by selection, many results point to other non-random factors that play a role in karyotype evolution. In cancer, chromosomal instability leads to characteristic changes in the karyotype, in which different individuals with a specific type of cancer display similar changes in karyotype structure over time. Statistical analyses of chromosome lengths in different species indicate that the length distribution of chromosomes is not consistent with models in which the lengths of chromosomes are random or evolve solely by simple random processes. A better understanding of the mechanisms underlying karyotype evolution should enable the development of quantitative theoretical models that combine the random and deterministic processes that can be compared to experimental determinations of the karyotype in diverse settings.

Topology and Dynamics of Transcriptome (Dys)Regulation.

RNA transcripts play a crucial role as witnesses of gene expression health. Identifying disruptive short sequences in RNA transcription and regulation is essential for potentially treating diseases. Let us delve into the mathematical intricacies of these sequences. We have previously devised a mathematical approach for defining a "healthy" sequence. This sequence is characterized by having at most four distinct nucleotides (denoted as nt≤4). It serves as the generator of a group denoted as fp. The desired properties of this sequence are as follows: fp should be close to a free group of rank nt-1, it must be aperiodic, and fp should not have isolated singularities within its SL2(C) character variety (specifically within the corresponding Groebner basis). Now, let us explore the concept of singularities. There are cubic surfaces associated with the character variety of a four-punctured sphere denoted as S24. When we encounter these singularities, we find ourselves dealing with some algebraic solutions of a dynamical second-order differential (and transcendental) equation known as the Painlevé VI Equation. In certain cases, S24 degenerates, in the sense that two punctures collapse, resulting in a "wild" dynamics governed by the Painlevé equations of an index lower than VI. In our paper, we provide examples of these fascinating mathematical structures within the context of miRNAs. Specifically, we find a clear relationship between decorated character varieties of Painlevé equations and the character variety calculated from the seed of oncomirs. These findings should find many applications including cancer research and the investigation of neurodegenative diseases.

SEEI: spherical evolution with feedback mechanism for identifying epistatic interactions.

BACKGROUND: Detecting epistatic interactions (EIs) involves the exploration of associations among single nucleotide polymorphisms (SNPs) and complex diseases, which is an important task in genome-wide association studies. The EI detection problem is dependent on epistasis models and corresponding optimization methods. Although various models and methods have been proposed to detect EIs, identifying EIs efficiently and accurately is still a challenge. RESULTS: Here, we propose a linear mixed statistical epistasis model (LMSE) and a spherical evolution approach with a feedback mechanism (named SEEI). The LMSE model expands the existing single epistasis models such as LR-Score, K2-Score, Mutual information, and Gini index. The SEEI includes an adaptive spherical search strategy and population updating strategy, which ensures that the algorithm is not easily trapped in local optima. We analyzed the performances of 8 random disease models, 12 disease models with marginal effects, 30 disease models without marginal effects, and 10 high-order disease models. The 60 simulated disease models and a real breast cancer dataset were used to evaluate eight algorithms (SEEI, EACO, EpiACO, FDHEIW, MP-HS-DHSI, NHSA-DHSC, SNPHarvester, CSE). Three evaluation criteria (pow1, pow2, pow3), a T-test, and a Friedman test were used to compare the performances of these algorithms. The results show that the SEEI algorithm (order 1, averages ranks = 13.125) outperformed the other algorithms in detecting EIs. CONCLUSIONS: Here, we propose an LMSE model and an evolutionary computing method (SEEI) to solve the optimization problem of the LMSE model. The proposed method performed better than the other seven algorithms tested in its ability to identify EIs in genome-wide association datasets. We identified new SNP-SNP combinations in the real breast cancer dataset and verified the results. Our findings provide new insights for the diagnosis and treatment of breast cancer. AVAILABILITY AND IMPLEMENTATION: https://github.com/scutdy/SSO/blob/master/SEEI.zip .

Synergistic epistasis among cancer drivers can rescue early tumors from the accumulation of deleterious passengers.

Epistasis among driver mutations is pervasive and explains relevant features of cancer, such as differential therapy response and convergence towards well-characterized molecular subtypes. Furthermore, a growing body of evidence suggests that tumor development could be hampered by the accumulation of slightly deleterious passenger mutations. In this work, we combined empirical epistasis networks, computer simulations, and mathematical models to explore how synergistic interactions among driver mutations affect cancer progression under the burden of slightly deleterious passengers. We found that epistasis plays a crucial role in tumor development by promoting the transformation of precancerous clones into rapidly growing tumors through a process that is analogous to evolutionary rescue. The triggering of epistasis-driven rescue is strongly dependent on the intensity of epistasis and could be a key rate-limiting step in many tumors, contributing to their unpredictability. As a result, central genes in cancer epistasis networks appear as key intervention targets for cancer therapy.

Mathematical model for the role of multiple pericentromeric repeats on heterochromatin assembly.

Although the length and constituting sequences for pericentromeric repeats are highly variable across eukaryotes, the presence of multiple pericentromeric repeats is one of the conserved features of the eukaryotic chromosomes. Pericentromeric heterochromatin is often misregulated in human diseases, with the expansion of pericentromeric repeats in human solid cancers. In this article, we have developed a mathematical model of the RNAi-dependent methylation of H3K9 in the pericentromeric region of fission yeast. Our model, which takes copy number as an explicit parameter, predicts that the pericentromere is silenced only if there are many copies of repeats. It becomes bistable or desilenced if the copy number of repeats is reduced. This suggests that the copy number of pericentromeric repeats alone can determine the fate of heterochromatin silencing in fission yeast. Through sensitivity analysis, we identified parameters that favor bistability and desilencing. Stochastic simulation shows that faster cell division and noise favor the desilenced state. These results show the unexpected role of pericentromeric repeat copy number in gene silencing and provide a quantitative basis for how the copy number allows or protects repetitive and unique parts of the genome from heterochromatin silencing, respectively.

A phylogenetic method linking nucleotide substitution rates to rates of continuous trait evolution.

Genomes contain conserved non-coding sequences that perform important biological functions, such as gene regulation. We present a phylogenetic method, PhyloAcc-C, that associates nucleotide substitution rates with changes in a continuous trait of interest. The method takes as input a multiple sequence alignment of conserved elements, continuous trait data observed in extant species, and a background phylogeny and substitution process. Gibbs sampling is used to assign rate categories (background, conserved, accelerated) to lineages and explore whether the assigned rate categories are associated with increases or decreases in the rate of trait evolution. We test our method using simulations and then illustrate its application using mammalian body size and lifespan data previously analyzed with respect to protein coding genes. Like other studies, we find processes such as tumor suppression, telomere maintenance, and p53 regulation to be related to changes in longevity and body size. In addition, we also find that skeletal genes, and developmental processes, such as sprouting angiogenesis, are relevant.

A novel hypergraph model for identifying and prioritizing personalized drivers in cancer.

Cancer development is driven by an accumulation of a small number of driver genetic mutations that confer the selective growth advantage to the cell, while most passenger mutations do not contribute to tumor progression. The identification of these driver genes responsible for tumorigenesis is a crucial step in designing effective cancer treatments. Although many computational methods have been developed with this purpose, the majority of existing methods solely provided a single driver gene list for the entire cohort of patients, ignoring the high heterogeneity of driver events across patients. It remains challenging to identify the personalized driver genes. Here, we propose a novel method (PDRWH), which aims to prioritize the mutated genes of a single patient based on their impact on the abnormal expression of downstream genes across a group of patients who share the co-mutation genes and similar gene expression profiles. The wide experimental results on 16 cancer datasets from TCGA showed that PDRWH excels in identifying known general driver genes and tumor-specific drivers. In the comparative testing across five cancer types, PDRWH outperformed existing individual-level methods as well as cohort-level methods. Our results also demonstrated that PDRWH could identify both common and rare drivers. The personalized driver profiles could improve tumor stratification, providing new insights into understanding tumor heterogeneity and taking a further step toward personalized treatment. We also validated one of our predicted novel personalized driver genes on tumor cell proliferation by vitro cell-based assays, the promoting effect of the high expression of Low-density lipoprotein receptor-related protein 1 (LRP1) on tumor cell proliferation.

polyGBLUP: a modified genomic best linear unbiased prediction improved the genomic prediction efficiency for autopolyploid species.

Given the universality of autopolyploid species in nature, it is crucial to develop genomic selection methods that consider different allele dosages for autopolyploid breeding. However, no method has been developed to deal with autopolyploid data regardless of the ploidy level. In this study, we developed a modified genomic best linear unbiased prediction (GBLUP) model (polyGBLUP) through constructing additive and dominant genomic relationship matrices based on different allele dosages. polyGBLUP could carry out genomic prediction for autopolyploid species regardless of the ploidy level. Through comprehensive simulations and analysis of real data of autotetraploid blueberry and guinea grass and autohexaploid sweet potato, the results showed that polyGBLUP achieved higher prediction accuracy than GBLUP and its superiority was more obvious when the ploidy level of autopolyploids is high. Furthermore, when the dominant effect was added to polyGBLUP (polyGDBLUP), the greater the dominance degree, the more obvious the advantages of polyGDBLUP over the diploid models in terms of prediction accuracy, bias, mean squared error and mean absolute error. For real data, the superiority of polyGBLUP over GBLUP appeared in blueberry and sweet potato populations and a part of the traits in guinea grass population due to the high correlation coefficients between diploid and polyploidy genomic relationship matrices. In addition, polyGDBLUP did not produce higher prediction accuracy than polyGBLUP for most traits of real data as dominant genetic variance was not captured for these traits. Our study will be a significant promising method for genomic prediction of autopolyploid species.

Distributions of 4-subtree patterns for uniform random unrooted phylogenetic trees.

Tree shape statistics based on peripheral structures have been utilized to study evolutionary mechanisms and inference methods. Partially motivated by a recent study by Pouryahya and Sankoff on modeling the accumulation of subgenomes in the evolution of polyploids, we present the distribution of subtree patterns with four or fewer leaves for the unrooted Proportional to Distinguishable Arrangements (PDA) model. We derive a recursive formula for computing the joint distributions, as well as a Strong Law of Large Numbers and a Central Limit Theorem for the joint distributions. This enables us to confirm several conjectures proposed by Pouryahya and Sankoff, as well as provide some theoretical insights into their observations. Based on their empirical datasets, we demonstrate that the statistical test based on the joint distribution could be more sensitive than those based on one individual subtree pattern to detect the existence of evolutionary forces such as whole genome duplication.

IGF2BP2 and IGFBP3 Genotypes, Haplotypes, and Genetic Models Studies in Polycystic Ovary Syndrome.

BACKGROUND: Insulin resistance has been correlated with the genetic diversity within the insulin-like binding proteins genes. Moreover, insulin resistance is one of the key characteristics of the widespread reproductive endocrine condition known as polycystic ovarian syndrome (PCOS). Hence, this study is aimed to determine the association between IGFBP3 and IGF2BP2 gene variants and PCOS risk. METHODS: A total of 300 subjects (150 PCOS cases diagnosed based on Rotterdam ESHRE/ASRM consensus criteria and 150 healthy subjects) were recruited in this case-control cross-sectional study. Tetra-primer amplification refractory mutation system polymerase chain reaction (ARMS-PCR) was used for genotyping rs11705701, whereas genotyping of rs1470579 and rs2854744 was done employing PCR-restriction fragment length polymorphism (PCR-RFLP) technique. RESULTS: The CC and AA+AC genotypes of rs1470579 conferred an increased risk of PCOS in our population. Regarding the rs2854744, an increased risk of PCOS was observed under the codominant homozygous (TT vs. GG) model by 2.54 fold. The C allele of rs1470579 and T allele of rs2854744 enhanced PCOS risk by 1.97 and 1.46 folds, respectively. Haplotype analysis showed that the Ars1470579Ars11705701 haplotype conferred a decreased risk of PCOS (odds ratio = 0.53, 95% confidence interval = 0.34-0.83, p = 0.006). The AC/GG/GT, AA/GA/GT, AC/GA/GG, and AC/GA/GT genotype combinations of rs1470579/rs11705701/rs2854744 were associated with a decreased risk of the disease. CONCLUSIONS: IGF2BP2 rs1470579 and IGFBP3 rs2854744 enhanced PCOS susceptibility in a Southeastern Iranian population. Further investigation involving larger cohorts representing diverse ethnic backgrounds is needed to confirm the current findings.

Revealing genomic heterogeneity and commonality: A penalized integrative analysis approach accounting for the adjacency structure of measurements.

Advancements in high-throughput genomic technologies have revolutionized the field of disease biomarker identification by providing large-scale genomic data. There is an increasing focus on understanding the relationships among diverse patient groups with distinct disease subtypes and characteristics. Complex diseases exhibit both heterogeneity and shared genomic factors, making it essential to investigate these patterns to accurately detect markers and comprehensively understand the diseases. Integrative analysis has emerged as a promising approach to address this challenge. However, existing studies have been limited by ignoring the adjacency structure of genomic measurements, such as single nucleotide polymorphisms (SNPs) and DNA methylations. In this study, we propose a structured integrative analysis method that incorporates a spline type penalty to accommodate this adjacency structure. We utilize a fused lasso type penalty to identify both heterogeneity and commonality across the groups. Extensive simulations demonstrate its superiority compared to several direct competing methods. The analysis of The Cancer Genome Atlas melanoma data with DNA methylation measurements and GENEVA diabetes data with SNP measurements exhibit that the proposed analysis lead to meaningful findings with better prediction performance and higher selection stability.

Optimizing therapeutic targets for breast cancer using boolean network models.

Studying gene regulatory networks associated with cancer provides valuable insights for therapeutic purposes, given that cancer is fundamentally a genetic disease. However, as the number of genes in the system increases, the complexity arising from the interconnections between network components grows exponentially. In this study, using Boolean logic to adjust the existing relationships between network components has facilitated simplifying the modeling process, enabling the generation of attractors that represent cell phenotypes based on breast cancer RNA-seq data. A key therapeutic objective is to guide cells, through targeted interventions, to transition from the current cancer attractor to a physiologically distinct attractor unrelated to cancer. To achieve this, we developed a computational method that identifies network nodes whose inhibition can facilitate the desired transition from one tumor attractor to another associated with apoptosis, leveraging transcriptomic data from cell lines. To validate the model, we utilized previously published in vitro experiments where the downregulation of specific proteins resulted in cell growth arrest and death of a breast cancer cell line. The method proposed in this manuscript combines diverse data sources, conducts structural network analysis, and incorporates relevant biological knowledge on apoptosis in cancer cells. This comprehensive approach aims to identify potential targets of significance for personalized medicine.

Joint Linkage and Association Analysis Using GENEHUNTER-MODSCORE with an Application to Familial Pancreatic Cancer.

INTRODUCTION: Joint linkage and association (JLA) analysis combines two disease gene mapping strategies: linkage information contained in families and association information contained in populations. Such a JLA analysis can increase mapping power, especially when the evidence for both linkage and association is low to moderate. Similarly, an association analysis based on haplotypes instead of single markers can increase mapping power when the association pattern is complex. METHODS: In this paper, we present an extension to the GENEHUNTER-MODSCORE software package that enables a JLA analysis based on haplotypes and uses information from arbitrary pedigree types and unrelated individuals. Our new JLA method is an extension of the MOD score approach for linkage analysis, which allows the estimation of trait-model and linkage disequilibrium (LD) parameters, i.e., penetrance, disease-allele frequency, and haplotype frequencies. LD is modeled between alleles at a single diallelic disease locus and up to three diallelic test markers. Linkage information is contributed by additional multi-allelic flanking markers. We investigated the statistical properties of our JLA implementation using extensive simulations, and we compared our approach to another commonly used single-marker JLA test. To demonstrate the applicability of our new method in practice, we analyzed pedigree data from the German National Case Collection for Familial Pancreatic Cancer (FaPaCa). RESULTS: Based on the simulated data, we demonstrated the validity of our JLA-MOD score analysis implementation and identified scenarios in which haplotype-based tests outperformed the single-marker test. The estimated trait-model and LD parameters were in good accordance with the simulated values. Our method outperformed another commonly used JLA single-marker test when the LD pattern was complex. The exploratory analysis of the FaPaCa families led to the identification of a promising genetic region on chromosome 22q13.33, which can serve as a starting point for future mutation analysis and molecular research in pancreatic cancer. CONCLUSION: Our newly proposed JLA-MOD score method proves to be a valuable gene mapping and characterization tool, especially when either linkage or association information alone provide insufficient power to identify the disease-causing genetic variants.

A novel genetic model provides a unique perspective on the relationship between postexercise glycogen concentration and increases in the abundance of key metabolic proteins after acute exercise.

Some acute exercise effects are influenced by postexercise (PEX) diet, and these diet-effects are attributed to differential glycogen resynthesis. However, this idea is challenging to test rigorously. Therefore, we devised a novel genetic model to modify muscle glycogen synthase 1 (GS1) expression in rat skeletal muscle with an adeno-associated virus (AAV) short hairpin RNA knockdown vector targeting GS1 (shRNA-GS1). Contralateral muscles were injected with scrambled shRNA (shRNA-Scr). Muscles from exercised (2-hour-swim) and time-matched sedentary (Sed) rats were collected immediately postexercise (IPEX), 5-hours-PEX (5hPEX), or 9-hours-PEX (9hPEX). Rats in 5hPEX and 9hPEX experiments were refed (RF) or not-refed (NRF) chow. Muscles were analyzed for glycogen, abundance of metabolic proteins (pyruvate dehydrogenase kinase 4, PDK4; peroxisome proliferator-activated receptor γ coactivator-1α, PGC1α; hexokinase II, HKII; glucose transporter 4, GLUT4), AMP-activated protein kinase phosphorylation (pAMPK), and glycogen metabolism-related enzymes (glycogen phosphorylase, PYGM; glycogen debranching enzyme, AGL; glycogen branching enzyme, GBE1). shRNA-GS1 versus paired shRNA-Scr muscles had markedly lower GS1 abundance. IPEX versus Sed rats had lower glycogen and greater pAMPK, and neither of these IPEX-values differed for shRNA-GS1 versus paired shRNA-Scr muscles. IPEX versus Sed groups did not differ for abundance of metabolic proteins, regardless of GS1 knockdown. Glycogen in RF-rats was lower for shRNA-GS1 versus paired shRNA-Scr muscles at both 5hPEX and 9hPEX. HKII protein abundance was greater for 5hPEX versus Sed groups, regardless of GS1 knockdown or diet, and despite differing glycogen levels. At 9hPEX, shRNA-GS1 versus paired shRNA-Scr muscles had greater PDK4 and PGC1α abundance within each diet group. However, the magnitude of PDK4 or PGC1α changes was similar in each diet group regardless of GS1 knockdown although glycogen differed between paired muscles only in RF-rats. In summary, we established a novel genetic approach to investigate the relationship between muscle glycogen and other exercise effects. Our results suggest that exercise-effects on abundance of several metabolic proteins did not uniformly correspond to differences in postexercise glycogen.

Haplotype reconstruction for genetically complex regions with ambiguous genotype calls: Illustration by the KIR gene region.

Advances in DNA sequencing technologies have enabled genotyping of complex genetic regions exhibiting copy number variation and high allelic diversity, yet it is impossible to derive exact genotypes in all cases, often resulting in ambiguous genotype calls, that is, partially missing data. An example of such a gene region is the killer-cell immunoglobulin-like receptor (KIR) genes. These genes are of special interest in the context of allogeneic hematopoietic stem cell transplantation. For such complex gene regions, current haplotype reconstruction methods are not feasible as they cannot cope with the complexity of the data. We present an expectation-maximization (EM)-algorithm to estimate haplotype frequencies (HTFs) which deals with the missing data components, and takes into account linkage disequilibrium (LD) between genes. To cope with the exponential increase in the number of haplotypes as genes are added, we add three components to a standard EM-algorithm implementation. First, reconstruction is performed iteratively, adding one gene at a time. Second, after each step, haplotypes with frequencies below a threshold are collapsed in a rare haplotype group. Third, the HTF of the rare haplotype group is profiled in subsequent iterations to improve estimates. A simulation study evaluates the effect of combining information of multiple genes on the estimates of these frequencies. We show that estimated HTFs are approximately unbiased. Our simulation study shows that the EM-algorithm is able to combine information from multiple genes when LD is high, whereas increased ambiguity levels increase bias. Linear regression models based on this EM, show that a large number of haplotypes can be problematic for unbiased effect size estimation and that models need to be sparse. In a real data analysis of KIR genotypes, we compare HTFs to those obtained in an independent study. Our new EM-algorithm-based method is the first to account for the full genetic architecture of complex gene regions, such as the KIR gene region. This algorithm can handle the numerous observed ambiguities, and allows for the collapsing of haplotypes to perform implicit dimension reduction. Combining information from multiple genes improves haplotype reconstruction.

BRASS: Permutation methods for binary traits in genetic association studies with structured samples.

In genetic association analysis of complex traits, permutation testing can be a valuable tool for assessing significance when the distribution of the test statistic is unknown or not well-approximated. This commonly arises, e.g, in tests of gene-set, pathway or genome-wide significance, or when the statistic is formed by machine learning or data adaptive methods. Existing applications include eQTL mapping, association testing with rare variants, inclusion of admixed individuals in genetic association analysis, and epistasis detection among many others. For genetic association testing in samples with population structure and/or relatedness, use of naive permutation can lead to inflated type 1 error. To address this in quantitative traits, the MVNpermute method was developed. However, for association mapping of a binary trait, the relationship between the mean and variance makes both naive permutation and the MVNpermute method invalid. We propose BRASS, a permutation method for binary traits, for use in association mapping in structured samples. In addition to modeling structure in the sample, BRASS allows for covariates, ascertainment and simultaneous testing of multiple markers, and it accommodates a wide range of test statistics. In simulation studies, we compare BRASS to other permutation and resampling-based methods in a range of scenarios that include population structure, familial relatedness, ascertainment and phenotype model misspecification. In these settings, we demonstrate the superior control of type 1 error by BRASS compared to the other 6 methods considered. We apply BRASS to assess genome-wide significance for association analyses in domestic dog for elbow dysplasia (ED) and idiopathic epilepsy (IE). For both traits we detect previously identified associations, and in addition, for ED, we detect significant association with a SNP on chromosome 35 that was not detected by previous analyses, demonstrating the potential of the method.

Identification of Genomic Regions Implicated in Susceptibility to Schistosoma mansoni Infection in a Murine Backcross Genetic Model.

Only a small number of infected people are highly susceptible to schistosomiasis, showing high levels of infection or severe liver fibrosis. The susceptibility to schistosome infection is influenced by genetic background. To assess the genetic basis of susceptibility and identify the chromosomal regions involved, a backcross strategy was employed to generate high variation in schistosomiasis susceptibility. This strategy involved crossing the resistant C57BL/6J mouse strain with the susceptible CBA/2J strain. The resulting F1 females (C57BL/6J × CBA/2J) were then backcrossed with CBA/2J males to generate the backcross (BX) cohort. The BX mice exhibited a range of phenotypes, with disease severity varying from mild to severe disease, lacking a fully resistant group. We observed four levels of infection intensity using cluster and principal component analyses and K-means based on parasitological, pathological, and immunological trait measurements. The mice were genotyped with 961 informative SNPs, leading to the identification of 19 new quantitative trait loci (QTL) associated with parasite burden, liver lesions, white blood cell populations, and antibody responses. Two QTLs located on chromosomes 15 and 18 were linked to the number of granulomas, liver lesions, and IgM levels. The corresponding syntenic human regions are located in chromosomes 8 and 18. None of the significant QTLs had been reported previously.

Bioactivity-Guided Isolation and Antihypertensive Activity of Citrullus colocynthis Polyphenols in Rats with Genetic Model of Hypertension.

Background and Objectives: Citrullus colocynthis belongs to the Cucurbitaceae family and is a wild medicinal plant used in folk literature to treat various diseases. The purpose of the current study was to explore the antihypertensive and antioxidant potentials of Citrullus colocynthis (CC) polyphenol-rich fractions using a spontaneous hypertensive rat (SHR) model. Materials and Methods: The concentrated aqueous ethanol extract of CC fruit was successively fractioned using solvents of increasing polarity, i.e., hexane, chloroform, ethyl acetate and n-butanol. The obtained extracts were analyzed for total phenolic content (TPC), total flavonoid content (TFC) and total flavonol content (TOF). Moreover, the CC extracts were further evaluated for radical scavenging capacity using 2,2-diphenyl-1-picrylhydrazyl (DPPH) and 2,2-azino-bis-3-ethylbenzothiazoline-6-sulfonic acid (ABTS) assays and antioxidant activity using inhibition of linoleic acid peroxidation and determination of reducing potential protocols. The phytochemical components were characterized by HPLC-MWD-ESI-MS in positive ionization mode. Results: The results showed that ethyl acetate fraction (EAF) exhibited a higher content of phenolic compounds in term of TPC (289 mg/g), TFC (7.6 mg/g) and TOF (35.7 mg/g). EAF showed higher antioxidant and DPPH and ABTS scavenging activities with SC50 values of 6.2 and 79.5 µg/mL, respectively. LCMS analysis revealed that twenty polyphenol compounds were identified in the EAF, including phenolic acids and flavonoids, mainly myricetin and quercetin derivatives. The in vivo antihypertensive activity of EAF of CC on SHR revealed that it significantly decreased the mean arterial pressure (MAP), systolic blood pressure (SBP), diastolic blood pressures (DBP) and pulse pressure (PP) as compared to normal and hypertensive control groups. Moreover, EAF of CC significantly reduced the oxidative stress in the animals in a dose-dependent manner by normalizing the levels of superoxide dismutase (SOD), malondialdehyde (MDA), reduced glutathione (GSH), nitric oxide (NOx) and total antioxidant capacity (TAC). Furthermore, the treatment groups, especially the 500 mg of EAF per kg body weight (EA-500) group, significantly (p ≤ 0.05) improved the electrocardiogram (ECG) pattern and pulse wave velocity (PWV). Conclusion: It was concluded that the EAF of CC is a rich source of polyphenols and showed the best antioxidant activity and antihypertensive potential in SHR.

Data-adaptive and pathway-based tests for association studies between somatic mutations and germline variations in human cancers.

Cancer is a disease driven by a combination of inherited genetic variants and somatic mutations. Recently available large-scale sequencing data of cancer genomes have provided an unprecedented opportunity to study the interactions between them. However, previous studies on this topic have been limited by simple, low statistical power tests such as Fisher's exact test. In this paper, we design data-adaptive and pathway-based tests based on the score statistic for association studies between somatic mutations and germline variations. Previous research has shown that two single-nucleotide polymorphism (SNP)-set-based association tests, adaptive sum of powered score (aSPU) and data-adaptive pathway-based (aSPUpath) tests, increase the power in genome-wide association studies (GWASs) with a single disease trait in a case-control study. We extend aSPU and aSPUpath to multi-traits, that is, somatic mutations of multiple genes in a cohort study, allowing extensive information aggregation at both SNP and gene levels. p $p$ -values from different parameters assuming varying genetic architecture are combined to yield data-adaptive tests for somatic mutations and germline variations. Extensive simulations show that, in comparison with some commonly used methods, our data-adaptive somatic mutations/germline variations tests can be applied to multiple germline SNPs/genes/pathways, and generally have much higher statistical powers while maintaining the appropriate type I error. The proposed tests are applied to a large-scale real-world International Cancer Genome Consortium whole genome sequencing data set of 2583 subjects, detecting more significant and biologically relevant associations compared with the other existing methods on both gene and pathway levels. Our study has systematically identified the associations between various germline variations and somatic mutations across different cancer types, which potentially provides valuable utility for cancer risk prediction, prognosis, and therapeutics.

ioSearch: An approach for identifying interacting multiomics biomarkers using a novel algorithm with application on breast cancer data sets.

Identification of biomarkers by integrating multiple omics together is important because complex diseases occur due to an intricate interplay of various genetic materials. Traditional single-omics association tests neither explore this crucial interomics dependence nor identify moderately weak signals due to the multiple-testing burden. Conversely, multiomics data integration imparts complementary information but suffers from an increased multiple-testing burden, data diversity inherent with different omics features, high-dimensionality, and so forth. Most of the available methods address subtype classification using dimension-reduction techniques to circumvent the sample size issue but interacting multiomics biomarker identification methods are unavailable. We propose a two-step model that first investigates phenotype-omics association using logistic regression. Then, selects disease-associated omics using sparse principal components which explores the interrelationship of multiple variables from two omics in a multivariate multiple regression framework. On the basis of this model, we developed a multiomics biomarker identification algorithm, interacting omics search (ioSearch), that jointly tests the effect of multiple omics with disease and between-omics associations by using pathway information that subsequently reduces the multiple-testing burden. Further, inference in terms of p values potentially makes it an easily interpretable biomarker identification tool. Extensive simulation demonstrates ioSearch as statistically powerful with a controlled Type-I error rate. Its application to publicly available breast cancer data sets identified relevant omics features in important pathways.